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Calculation of quantum tunneling for a spatially extended defect: the dislocation kink 

in copper has a low effective mass. 
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tLaboratory of Atomic and Solid State Physics (LASSP), Clark Hall,Cornell University, Ithaca, NY 14853-2501, USA 
^Cornell Theory Center, Cornell University, Ithaca, NY 14853, USA 

Several experiments indicate that there are atomic tunneling defects in plastically deformed metals. 
How this is possible has not been clear, given the large mass of the metal atoms. Using a classical 
molecular-dynamics calculation, we determine the structures, energy barriers, effective masses, and 
quantum tunneling rates for dislocation kinks and jogs in copper screw dislocations. We find that jogs 
are unlikely to tunnel, but the kinks should have large quantum fluctuations. The kink motion involves 
hundreds of atoms each shifting a tiny amount, leading to a small effective mass and tunneling barrier. 
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Tunneling of atoms is unusual. At root, the reason 
atoms don't tunnel is that their tunneling barriers and 
distances are set by the much lighter electrons. The tun- 
neling of a proton over a barrier one Rydberg high and 
one Bohr radius wide is suppressed by the exponential of 
y/2M p Ryal = ^/M p /m e ~ 42.85: a factor of 10~ 19 . 

Nonetheless, atomic quantum tunneling dominates the 
low temperature properties of glasses and many doped 
crystals Q. In glasses, there are rare regions (one per 
10 5 or 10 6 molecular units) where an atom or group of 
atoms has a double well with low enough barrier, tun- 
neling distance, and asymmetry to be active. For cer- 
tain dopants in crystals, off-center atoms and rotational 
modes of nearly spherical ionic molecules have unusu- 
ally low barriers and tunneling distances. Quantitative 
modeling of these spatially localized tunneling defects has 
been frustrated by the demands for extremely accurate 
estimates of energy barriers, beyond even the best den- 
sity functional electronic structure calculations available 
today. Also, although all detailed models of tunneling in 
glasses have basically involved one or a very few atoms, 
there has long been speculation that large numbers of 
atoms may be shifting during the tunneling process. || 

There is much evidence that quantum tunneling is 
important to the properties of undoped, plastically 
deformed metals. Quantum creep glassy low- 

temperature behavior || , and two-channel Kondo scaling 
seen in the voltage and temperature-dependent electri- 
cal conductivity in nanoconstrictions || have been at- 
tributed to quantum tunneling associated with disloca- 
tions. It has never been clear how this can occur, given 
the large masses of the metal atoms involved. 

We show here using a classical effective-medium inter- 
atomic potential that quantum fluctuations can indeed 
be important in the dynamics of one particular defect: a 
kink in the split-core screw dislocation in copper. The 
motion of the kink involves a concerted motion of hun- 
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FIG. 1. Kink geometry. Broadly speaking, the screw 
dislocation represents the locus where planes of atoms form a 
helix. In copper it spreads out into a ribbon along the x-axis 
to lower its energy. The kink we study shifts the ribbon by 
one atom in the x direction. More specifically, the b=| [110] 
dislocation on the left dissociates on the (111) plane into the 
Shockley partials: aC = f [121] and Dq = |[2lT] respectively. 
The kink is introduced with line vector Ckink ~ § [112], and 
dissociates into a wide screw-like- and a bulky edge-like kink 
located on the partial dislocations. The lighter atoms are on 
the stacking fault (hep local environments) and the darker 
atoms are along the partial dislocations (neither hep nor fee). 

dreds of copper atoms, leading to a dramatic decrease 
in its effective mass. This derealization perhaps lends 
support to ideas about collective centers in glasses. Also, 
because our important conclusions rest upon this dereal- 
ization they are qualitatively much less sensitive to the 
accuracy of our potential than calculations for spatially 
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localized tunneling defects. We assert that these kinks 
are likely the only candidate for quantum tunneling in 
pure fee metals. 

The kink simulation consists of two screw dislocations 
with opposite Burgers vectors b=±| [110], allowing pe- 
riodic boundary conditions giving us the perfect transla- 
tional invariance necessary to measure energy differences 
to the accuracy we need. The two dislocations are placed 
in different (111) planes separated by 20 (111) planes (4.4 
ran), see figure |l|. The system is 86 planes wide (19.3 
nm) in the two (non-orthogonal) directions, and extends 
44.5 b (11.4 nm) along the dislocations. We introduce 
kinks or jogs on the dislocations by applying skew peri- 
odic boundary conditions to system, i.e. we introduce a 
small mismatch in the dislocation cores at the interface 
to the next cell. The procedure also introduces a row of 
interstitial atoms between the kinked dislocations, which 
is subsequently removed from the system, leaving us with 
a total of 329 102 atoms. The kinks have a net line vector 
of £kink = |[Tl2], with a the lattice constant. 

To show how unusual the properties of the kink are, 
we also study the properties of a dislocation jog. The 
jog simulation, and the associated energy barrier calcu- 
lation, is similar and is described elsewhere 0. The el- 
ementary jog we study is introduced with a line vector 
oriented in the (111) glide plane of the screw disloca- 
tion, £j g = | [112], which then transforms into an obtuse 
lower energy configuration: f [112] -> f [101] + f [TTO]. 
This jog is expected to be the most mobile of the jogs, 
second only to the kink in mobility. 

We introduce the two kinked dislocations directly as 
Shockley partial dislocations, see figure |l|, and relax using 
the MD-min algorithm ||, using Effective Medium The- 
ory (EMT): a many-body classical potential which 
is computationally almost as fast as a pair potential, 
while still describing the elastic properties well. The 
elastic constants of the potential are: C\\ = 176.2 GPa, 
C u = 116.0 GPa and C 44 = 90.6 GPa with a Voigt av- 
erage shear modulus of /i = 66 GPa, and an intrinsic 
stacking fault energy of ^ = 31 mJ/m 2 . 

We present three quantities for the kink and jog: the 
Peierls-like barrier for migration along the dislocation, 
the effective mass, and an upper bound for the WKB 
factor suppressing quantum tunneling through that bar- 
rier. Since the motion of these defects involves several 
atoms moving in a coordinated fashion, we use instan- 
tons: the appropriate generalization of WKB analysis to 
many-dimensional configuration spaces An up- 

per bound for the tunneling matrix element is given by 
the effective mass approximation pT]Jl3| , 



M*(Q) = *^Mi(dq i /dQ)' 2 



(2) 



A < Hloq exp 



y/2M*(Q)V(Q)dQ/h 



(1) 



where loq is an attempt frequency, V*(Q) is the energy 
of the defect at position Q with the neighbors in their 
relaxed, minimum energy positions qi(Q), and 



is the effective mass of the defect incorporating the ki- 
netic energy of the surrounding atoms as they respond 
adiabatically to its motion. The effective mass approx- 
imation is usually excellent for atomic tunneling. The 
method is variational, so equations [I] and ^| remain up- 
per bounds using other assumptions about the tunneling 
path qi{Q) (such as the straight- line path between the 
two minima described below for the kinks). 

The difficulty of finding models for atomic tunneling 
is illustrated rather well by the properties of the jog we 
study. The barrier for migration was determined to be 
15 meV 0: lower than other jogs, or even than sur- 
face diffusion barriers calculated with the same poten- 
tial. The effective mass for the jog, estimated by sum- 
ming the squared displacement of the 200 atoms with 
largest motion, is Mjog — 0.36 Mq u : the jog is spa- 
tially localized (it doesn't disassociate into partials), with 
a few atoms in the core of the jog carrying most of 
the motion. The WKB tunneling matrix element for 
the jog to tunnel a distance Q = 2.5 A over a bar- 
rier V — 0.015 eV is suppressed by a factor of roughly 
cxp(-y / 2Afj* og t/Q/n) ~ 10~ 14 . Jogs don't tunnel much. 
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FIG. 2. The magnitude of the atomic displacement field 
as the two partial kinks move, along the cores of the two 
partial dislocations, as a function of the z-coordinate. The 
data points are fitted to 'periodic' Lorentzian distributions. 
The FWHM value of the two partial kinks are r e( jg e = 13b 
(bullets) and Tscrew = 21b (triangles). The center of masses 
of the two partial kinks are separated: the screw-like partial 
kink is shifted by Sz — —1.5b. 

For the kinks, we take a relaxed initial configuration 
and define a final configuration with the kink migrated 
by one lattice spacing along the dislocation. The final 
position for each atom is given by the position of the 
neighboring atom closest to the current position minus 
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the kink migration vector l m [g r = § [HO] which repre- 
sents the net motion of the kink. This automatically gives 
the correct relaxed final position, which is otherwise dif- 
ficult to locate given the extremely small barriers. The 
width of the kinks is the traditional name for their extent 
along the axis of the screw dislocation. We can measure 
this width by looking at the net displacement of atoms 
between the initial and final configurations. We find that 
the displacement field is localized into two partial kinks, 
localized on the partial dislocation cores, see figure |^. 
These two partial kinks are quite wide (FWHM of 13 b 
and 21 b). They differ because the partials are of mixed 
edge and screw character; it is known that the kink 
which forces a mixed dislocation towards the screw di- 
rection will be wider and have higher energy. This is 
wider than the w < 10b predicted for slip dislocations in 
closed-packed materials by Hirth and Lothe, and Seeger 
and Schiller using line tension models fl5|| . 

Notice that the maximum net distance moved by an 
atom during the kink motion in figure ^|is around 0.01 A. 
Summing the squares of all the atomic motions, and using 
equation |[ we find an effective mass ^ — M^/ldO 
within the straight-line path approximation. This re- 
markably small mass can be attributed to three factors. 
(1) The mass is decreased because the screw dislocation 
is split into two partial dislocations |ll|. (2) The cores of 
the partial dislocations are spread transversally among 
Wt — 4<i, figure pj; this factor seems to have been missed 
in continuum treatments. These first two factors each 
reduce the total distance moved by an atom as the kink 
passes from z — —oo to +oo. (3) The kink partials aver- 
age Wt, — 17b wide (above), so the total atomic motion 
is spread between around 17 kink migration hops [ p"7j ]. 
Thus when the kink moves by x, the atoms in two regions 
I/Wl long and 1/Wt wide each move by x/(2W~lW~t), 
reducing [Tj]] the effective mass by roughly a factor of 
2W L W T ~ 136. 

Evaluating the energy at equally spaced atomic con- 
figurations and linearly interpolating between the initial 
and final states (along the straight-line path) yields an 
upper bound to the kink-migration barrier of 0.15 fJ-eV, 
figure |j. We attribute this extremely small barrier to 
the wide kink partials: we expect the barrier V to scale 
exponentially with the ratio of the kink width W to the 
(110) interplanar distance b: V ~ exp(— VT/b). If one 
thinks of the contribution to the energy of the n th layer 
as some analytic function f(n/W + 5), then the barrier 
is given by the variation of the sum J2^L-x> f( n /W + ^) 
with the position shift 6. The difference between this 
sum and the ((^-independent) integral is easily estimated 
by Fourier transforms, and is approximately 2f(2nW/a). 
The Fourier transform of an analytic function decays ex- 
ponentially. One imagines this could be proven to all 
orders in perturbation theory. 

This small barrier is not only negligible for thermal 
activation (two mK), but also for quantum tunneling. 



The WKB factor suppressing the tunneling would be 
exp{-^2M* hik VQ/h) = exp(-0.0148) = 0.985. Even 

at zero temperature, the kinks effectively act as free par- 
ticles, as suggested in the literature ( Q, among others). 
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FIG. 3. The magnitude of the atomic displacement field 
as the two partial kinks move, along the cores of the two 
partial dislocations, as a function of the x-coordinate. The 
core regions are fitted to two squared Lorentzian distributions. 
The partial core widths Wt — 4d (4 (111) lattice planes), 
significantly reduces the effective mass of the kink. 

Our estimated kink migration barrier is thus 10 5 times 
smaller than that for the most mobile of the jogs. How 
much can we trust our calculation of this remarkably 
small barrier? Schottky jll| estimates using a simple 
line-tension model that the barrier would be ~ 3 10 -5 eV 
in fee materials, using a Peierls stress er P = 10~ 2 /i and a 
kink width w — 10 b. This value is a factor of 200 higher 
than the barrier we find. On the other hand, both exper- 
iments and theoretical estimates predict tr P ~ 5 10 _6 /i 
for Cu fl9f| , yielding barriers orders of magnitude lower 
than ours. The interatomic potentials we use do not 
take into account directional bonding. This is usually 
a good approximation for noble metals; however, small 
contributions from angular forces may change the kink 
width. The kink width is like an energy barrier, balanc- 
ing different competing energies against one another: in 
analogy, we expect it to be accurate to within twenty or 
thirty percent. Our small value for the effective mass, 
dependent on the inverse cube of the spatial extent of 
the kink, is probably correct within a factor of two. The 
energy barrier is much more sensitive: if we take the to- 
tal exponential suppression to be 10 5 (using the jog as a 
"zero-length" defect) then each 20% change in the width 
would yield a factor of 10 change in the barrier height. 
The qualitative result of our calculation, that the barri- 
ers and effective masses are small, are robust not only to 
the use of an approximate classical potential, but may 
also apply to other noble metals and perhaps simple and 
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late transition metals. 
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FIG. 4. The activation energy as a function of the 
straight-line distance moved by the kink on one dislocation, 
with an associated barrier of -E ac t = 0.15 /ieV. Notice that 
this activation energy is about one part in 10 13 of the total 
system energy. 

In summary, we have used an atomistic calculation 
with classical potentials to extract energy barriers and 
effective masses for the quantum tunneling of dislocation 
jogs and kinks in copper. For jogs, the atomic displace- 
ments during tunneling arc primarily localized to a few 
atoms near the jog core, each moving a significant frac- 
tion of a lattice spacing. Consequently, the tunneling 
barrier and effective mass are relatively large, and tun- 
neling is unlikely. However, the kinks in screw disloca- 
tions are much more extended: as a kink moves by one 
lattice spacing, hundreds of atoms shift their positions 
by less than 1% of a lattice spacing. Both the energy 
barrier and the effective mass are reduced, to the extent 
that tunneling should occur readily. Kinks are likely the 
only candidate for quantum tunneling in pure crystalline 
materials. They may explain measurements of quantum 
creep, glassy internal friction, and non-magnetic Kondo 
effects seen in plastically deformed metals. 
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